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Abstract. Dimtri Mihalas' textbooks in the 70's and 80's on Stellar Atmospheres and Foundations 
of Radiation Hydrodynamics helped lay the early groundwork for understanding the moving atmo- 
spheres and winds of massive, luminous stars. Indeed, the central role of the momentum of stellar 
radiation in driving the mass outflow makes such massive-star winds key prototypes for radiation 
hydrodynamical processes. This paper reviews the dynamics of such radiative driving, building first 
upon the standard CAK model, and then discussing subtleties associated with the development and 
saturation of instabilities, and wind initiation near the sonic point base. An overall goal is to illumi- 
nate the rich physics of radiative driving and the challenges that lie ahead in developing dynamical 
models that can explain the broad scalings of mass loss rate and flow speed with stellar properties, 
as well as the often complex structure and variability observed in massive-star outflows. 
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INTRODUCTION 

Through his textbooks and numerous papers, Dimitri Mihalas helped lay the groundwork 
for radiation hydrodynamics as an active discipline within astrophysics. Many of the 
contributions to this conference in his honor have emphasized the common roles that 
radiation plays both as a diagnostic and as an energy source/sink for astrophysical fluids. 
But perhaps the most fundamental aspect of radiation hydrodynamics arises when the 
momentum of radiation plays the central role in the dynamical driving of an astrophysical 
flow. A key prototype for this lies in the strong stellar winds from hot, massive, luminous 
stars, which are driven by the scattering of the star's continuum radiation flux by line- 
transitions of metal ions. 

The discussion here aims to summarize some basic physical concepts from the nearly 
4 decades of research since line-driving was first proposed as the mechanism for hot- 
star winds [15]. I begin with a summary of the standard CAK/Sobolev [3, 38] formalism 
for spherical, steady-state models of such winds. I then review subsequent studies that 
emphasize how relaxation of the Sobolev approximation for localized line-transport is 
central to modeling both the strong, small-scale instability of line-driving, as well as 
the wind initiation near the transonic wind base. A subtle but important aspect of this is 
the strong, dynamical role of the diffuse component of line-scattered radiation, which is 
entirely ignored in the standard CAK/Sobolev approach, but which can be approximately 
accounted for through nonlocal, integral forms for the escape probability. 
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FIGURE 1. Two perspectives for the Doppler-shifted line-resonance in an accelerating flow. Right: 
Photons with a wavelength just shortward of a line propagate freely from the stellar surface up to a 
layer where the wind outflow Doppler shifts the line into a resonance over a narrow width (represented 
here by the shading) equal to the Sobolev length, set by the ratio of thermal speed to velocity gradient, 
= Vth/ {dv / dr). Left: Seen from successively larger radii within the accelerating wind, the Doppler- 
shift sweeps out an increasingly broadened line absorption trough in the stellar spectrum. 



THE CAK/SOBOLEV MODEL FOR STEADY WINDS 

Consider a steady-state stellar wind outflow in which the radiative acceleration gj-ad 
overcomes the local gravity GM^/r^ at radius r to drive a net acceleration v{dv/dr) in 
the radial flow speed v(r). Since overcoming gravity is key, it is convenient to define a 
dimensionless equation of motion that scales all accelerations by gravity, 

{l-Ws/w)w = -l+T,ad, (1) 

where Trad = grad / GM^, w = v^/v^^^., and w' = dw/dx, with x= l—R^/r and Vesc = 
^IGM^/R^ the escape speed from the stellar surface radius R^. Eqn. ([U) neglects gas 
pressure terms on the right side, since for isothermal sound speed a these are of order 
Ws = (a/vesc)^ ~ 0.001 compared to competing terms needed to drive the wind. 

For pure electron scattering opacity Kg, the scaled radiative acceleration is just the 
usual Eddington parameter 

r^^^^f^ = 2x10-5^2^. (2) 

Because typically ~ M^^^-^, stars with > IOMq have Eg > 10"^, with the Ed- 
dington limit Eg — > 1 perhaps even being central to setting a stellar upper mass limit 
of < ISOMq [9]. Eruptive mass loss from luminous blue variable (LBV) stars like 
T] Carinae might in fact be continuum-driven during episodes of super-Eddington lumi- 
nosity [|4|,|22l|28ll. 



But the resonant nature of line (bound-bound) scattering from metal ions leads to an 
opacity that is inherently much stronger than from free electrons. For example, in the 



somewhat idealized, optically thin limit that all the line opacity could be illuminated 
with a flat, unattenuated continuum spectrum with the full stellar luminosity, the total 
line-force would exceed the free-electron force by a huge factor, of order Q ^ 2000 



This implies line-driven winds can be initiated in even moderately massive stars 
with Fe > 5 X 10^^, while for more massive stars with Tg ^ I /2, the net outward line 
acceleration in principle could be as high as Tthin ~ Q^e ~ 1000 times the acceleration 
of gravity! 

In practice, self-absorption within strong lines limits the acceleration, with the mass 
loss rate M set at the level for which the line driving is just sufficient to overcome 
gravity. Indeed line- saturation keeps the dense, nearly static layers of the atmosphere 
gravitationally bound. But as illustrated by figure [H within the accelerating wind, the 
Doppler shift of the line-resonance out of the absorption shadow of underlying material 
exposes the line opacity to a less attenuated flux. This effectively desaturates the lines 
by limiting the resonance to a layer with width set by the Sobolev length, Isob — 
Vth/ [dv/dr), and with optical depth proportional iot= Kepc/ {dv/dr) = YgMc^ /L^w' . 

For the CAK line-ensemble with a power-law number distribution in line-strength, 
the cumulative force is reduced by a factor 1 / {Qt)^ from the optically thin value, 

TcAK = 7 , = F.ytr « = C(w') « , (3) 

(l-a)(eO" 

where the second equality defines the CAK "force muliplier" withQfc = "/( 1 — 
a). The last equality relates the line-force to the flow acceleration, with 
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Note that, for fixed sets of parameters for the star (L*, M*, Fg) and line-opacity (a, Q), 
this constant scales with the mass loss rate as C 1 /M". 

Neglecting the small sound-speed term Ws ~ 0.001 ^ 1, application of eqn. (O into 
(HI) gives the CAK equation of motion, 

F = w'-M-F,-C(w')'' = 0. (5) 

For small M (large C), there are two solutions, while for large M (small C), there are no 
solutions. The CAK critical solution corresponds to a maximal mass loss rate, defined 
by dF /dw' = 0, for which the C(w')" is tangent to the line 1 — Fg + at a critical 
acceleration = (1 — Te)a/{ \ — a). Since the scaled equation of motion ^ has no 
explicit spatial dependence, this critical acceleration applies throughout the wind, and 



Here we use a slight variation of the standard CAK notation in which the artificial dependence on 
a fiducial ion thermal speed is avoided by simply setting v,/, = c. Backconversion to CAK notation is 
achieved by multiplying t by v,/,/c and k by {v,i,/c)". The line normalization Q offers the advantages of 
being a dimensionless measure of line-opacity that is independent of the assumed ion thermal speed, with 
a nearly constant characteristic value of order 10^ for a wide range of ionization conditions 1 1 1 ] . 



so can be trivially integrated to yield w{x) = w[.x. In terms of dimensional quantities, 
this represents a specific case of the general "beta"-velocity-law, 
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v(r)=v.(^l-:^J , (6) 

where here /3 = 1/2, and the wind terminal speed Voo = Vesc^/ Oc{l — re)/(l — a). Sim- 
ilarly, the critical value Q yields, through eqn. (U), the standard CAK scaling for the 
mass loss rate 



(7) 



These CAK results strictly apply only under the idealized assumption that the stellar 
radiation is radially streaming from a point-source. If one takes into account the finite 
angular extent of the stellar disk, then near the stellar surface the radiative force is 
reduced by a factor /^/^ ^ 1/(1 + a), leading to a reduced mass loss rate [ilQ.i3Q1 

Mf, = f'^iTMcAK = ^Y^f^ - Mcak/2 . (8) 

Away from the star, the correction factor increases back toward unity, which for the re- 
duced base mass flux implies a stronger, more extended acceleration, giving a somewhat 
higher terminal speed, Voo ~ 3vesc, and a flatter velocity law, approximated by replacing 
the exponent in eqn. ©by /3 ^0.8. 

The effect of a radial change in ionization can be approximately taken into account by 
correcting the CAK force ([3]) by a factor of the form (ne/W) , where is the electron 
density, = 0.5(1 — -^Z 1 —R^/r) is the radiation "dilution factor", and the exponent has 
a typical value 5 ~ 0.1 lHH]. This factor introduces an additional density dependence to 
that already implied by the optical depth factor 1/?" given in eqn. ([3]). Its overall effect 
can be roughly accounted with the simple substition a ^ a' = a — 5 m the power 
exponents of the CAK mass loss scaling law CZ])- The general tendency is to moderately 
increase M, and accordingly to somewhat decrease the wind speed. 

The above scalings also ignore the finite gas pressure associated with a small but non- 
zero sound-speed parameter Ws- Through a perturbation expansion of the equation of 
motion ([T]) in this small parameter, it possible to derive simple scalings for the fractional 
corrections to the mass loss rate and terminal speed [27], 
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a Vesc 2(1 -a) ^l-aVesc 

For a typical case with a ^ 2/3 and w., = 0.001, the net effect is to increase the mass 
loss rate and decrease the wind terminal speed, both by about 10%. 

Finally, an important success of these CAK scaling laws is the theoretical rationale 
they provide for an empirically observed "Wind-Momentum-Luminosity" (WML) rela- 
tion rtl2[]. Combining the CAK mass-loss law ([7]) together with the scaling of the terminal 
speed with the effective escape, we obtain a WML relation of the form, 

Mv.cv^~Li/«'e'^"'~' (10) 



wherein we have neglected a residual dependence on M(l — Eg) that is generally very 
weak for the usual c ase that a' is near 2/3. Note that the direct dependence 2 ~ Z 
provides the scaling of the WML with metalicity Z. 



NON-SOBOLEV WIND INSTABILITY & INITIATION 

The above CAK steady-state model depends crucially on the use of the Sobolev approx- 
imation to compute the local CAK line force dB]). Analyses that relax this approximation 
show that the flow is subject to a strong, "line-deshadowing instability" (LDI) for ve- 
locity perturbations on a scale near and below the Sobolev length Isob = Vth/{dv/dr) 
12311 ■ Moreover, the diffuse, scattered component of the line force, which 



in the Sobolev limit is nullified by the fore-aft symmetry of the Sobolev escape proba- 
bility (see figure [lb), turns out to have important dynamics effects on both the instability 
and sonic -point initiation of the wind. 



Linear Analysis of Line-Deshadowing Instability 

For sinusoidal perturbations (~ with wavenumber k and frequency CO, the 

linearized momentum equation (ignoring the small gas pressure) relating the perturba- 
tions in velocity and radiative acceleration implies 

co = i^, (11) 
ov 



which shows that unstable growth, with 3(U > 0, requires 9t(5g/5v) > 0. For a purely 
Sobolev model [T], the CAK scaling of the line-force Q with velocity gradient v' 
implies 5g ~ 5v' ~ ikdv, giving a purely real CO, and thus a stable wave that propagates 
inward at phase speed, 

k dv 

which is now known as the "Abbott speed". Abbott [T] showed this is comparable to the 
outward wind flow speed, and in fact exactly equals it at the CAK critical point. 

As illustrated in figure [2^, instability arises from the deshadowing of the line by the 
extra Doppler shift from the velocity perturbation, giving 5g ~ dv and thus 3 ft) > 0. A 



general analysis [|25l1 yields a "bridging law" encompassing both effects. 



ov l+ikA 

where Q ~ gcak/vth sets the instability growth rate, and the "bridging length" A is found 
to be of order the Sobolev length l^ob- In the long-wavelength limit ^A ^ 1, we recover 
the stable, Abbott-wave scalings of the Sobolev approximation, dg/5v ^ ikQ.A = ikU; 
while in the short-wavelength limit fcA ^ 1, we obtain the instability scaling 5g ~ Q.dv. 
The instability growth rate is very large, about the flow rate through the Sobolev length. 
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FIGURE 2. (a) The line profile and direct intensity plotted vs. comoving frame frequency x — u = 
x — v/vfi,, with the light shaded overlap area proportional to the net direct line-force g^/,>. The dashed 
profile shows the effect of the Doppler shift from a perturbed velocity 5v, with the resulting extra area in 
the overlap with the blue-edge intensity giving a perturbed line-force 5g that scales in proportion to this 
perturbed velocity 8u — 5v/v(/,. (b) The comoving-frequency variation of the forward (+) and backward 
(-) streaming parts of the diffuse, scattered radiation, illustrating the cancelling of the overlap with the 
line profile that causes the net diffuse force to nearly vanish in a smooth, supersonic outflow, (c) However, 
because of the Doppler shift from the perturbed velocity, the dashed profile now has a stronger interaction 
with the backward streaming diffuse radiation. This results in a diffuse-line-drag force that scales with the 
negative of the perturbed velocity, and so tends to counter the instability of the dii-ect line-force in part a. 



Q ^ v/lsob- Since this is a large factor v/v//, bigger than the typical wind expansion rate 
dv / dr V / R^, a small perturbation at the wind base would, within this lineary theory, 
be amplified by an enormous factor, of order e^'/^r/, ~ giooi 



Numerical Simulations of Instability-Generated Wind Structure 



Numerical simulations of the nonlinear evolution require a non-Sobolev line-force 
computation on a spatial grid that spans the full wind expansion over several i?*, yet 
resolves the unstable structure at small scales near and below the Sobolev length. The 
first tractable approach [20] focussed on the absorption of the direct radiation from the 
stellar core, accounting now for the attenuation from intervening material by carrying 
out a nonlocal integral for the frequency-dependent radial optical depth, 

t{x,r) = r dr'K,p{/)(t> [x-vir')/vth] , (14) 

where ^ is the line-profile function, and x is the observer-frame frequency from line- 
center in units of the line thermal width. The corresponding nonlocal form for the CAK 
line-ensemble force from this direct stellar radiation is 

r.„.(r)^r.5'-°r''^ ^'^7 (") 

In the Sobolev approximation, t{x, r) ^ 4>(x — v/vth)t (where 4>(x) = j'^ ^{x!) dx'), this 
recovers the CAK form ([3]). But for perturbations on a spatial scale near and below 
the Sobolev length, its variation also scales in proportion to the perturbed velocity, 
leading to unstable amplification. Simulations show that because of inward nature of 
wave propagation implies an anti-correlation between velocity and density variation, the 
nonlinear growth leads to high-speed rarefactions that steepen into strong reverse shocks 
and compress material into dense clumps (or shells in these ID models) ^7^. 

The assumption of pure-absorption was criticized by Lucy [13], who pointed out 
that the interaction of a velocity perturbation with the background, diffuse radiation 
from line- scattering results in a line-drag effect that reduces, and potentially could even 
eliminate, the instability associated with the direct radiation from the underlying star. 
The basic effect is illustrated in figure [2l Panel b shows how the fore-aft (±) symmetry 
of the diffuse radiation leads to cancellation of the g+ and g force components from the 
forward and backward streams, as computed from a line-profile with frequency centered 
on the local comoving mean flow. In contrast, panel c shows that the Doppler shift 
associated with the velocity perturbation 5v breaks this symmetry, and leads to stronger 
forces from the component opposing the perturbation. 

Full linear stability analyses accounting for scattering effects 1 2^1231] show the frac- 



tion of the direct instability that is canceled by the line-drag of the perturbed diffuse 
force depends on the ratio of the scattering source function S to core intensity Ic, 

where the latter approximation applies for the optically thin form 2S/Ic = 1 — /i*. The 
net instability growth rate thus becomes 



This vanishes near the stellar surface, where /i* = 0, but it approaches half the pure- 
absorption rate far from the star, where /i* ^ 1 . This implies that the outer wind is still 
very unstable, with cumulative growth of ca. v^/2vth ~ 50 e-folds. 

Most efforts to account for scattering line-drag in simulations of the nonlinear evo- 
lution of the instability have centered on a Smooth Source Function (SSF) approach 
[1810,1231124']. This assumes that averaging over frequency and angle makes the scatter- 



ing source function relatively insensitive to flow structure, implying it can be pulled out 
of the integral in the formal solution for the diffuse intensity. Within a simple two-stream 
treatment of the line-transport, the net diffuse line-force then depends on the difference 
in the nonlocal escape probabilities b± associated with forward (-I-) vs. backward (-) in- 
tegrals of the frequency-dependent line-optical-depth (O. For a CAK line-ensemble, 
the net diffuse force can be written in a form quite analogous to the direct component 
([H, 

with 

whereH for f_ the integral bounds in (fT4l) are now from r to the outer radius R^ax ^M, [3] , 



and the overall normalization for T^iff assumes the optically thin source function from 
eqn. (fT6l) . In the Sobolev approximation, both the forward and backward integrals give 
the same form, viz. t±{±x,r) ^ (i>[±{x — v/vth)]t, leading to the net cancellation of the 
Sobolev diffuse force (figure[2b). But for perturbations on a spatial scale near and below 
the Sobolev length, the perturbed velocity breaks the forward/back symmetry (figure 
[2t), leading to perturbed diffuse force that now scales in proportion to the negative of the 
perturbed velocity, and thus giving the diffuse line-drag that reduces the net instability 
by the factors given in (fT6l) and (fTTl) . 

The left panel of figure [3] illustrates the results of a ID SSF simulation, starting from 
an initial condition set by smooth, steady-state CAK/Sobolev model (dashed curves). 
Because of the line-drag stabilization of the driving near the star (eqn.flTl). the wind base 
remains smooth and steady. But away from the stellar surface, the net strong instability 
leads to extensive structure in both velocity and density, roughly straddling the CAK 
steady-state. Because of the backstreaming component of the diffuse line-force causes 
any outer wind structure to induce small- amplitude fluctuations near the wind base, 
the wind structure, once initiated, is "self-excited", arising spontaneously without any 
explict perturbation from the stellar boundary. 

In the outer wind, the velocity variations become highly nonlinear and nonmonotonic, 
with amplitudes approaching 1000 km/s, leading to formation of strong shocks. How- 
ever, these high-velocity regions have very low density, and thus represent only very little 



To account for the nonradial radiation from a finite-angle star, in practice SSF models carry out the 
forward/back optical depth integrals along a ray with an impact parameter that is R^^j^T. off the disk 
center. A similar application for the direct force ( fTSl ) gives within a few percent the proper finite-disk 
correction to the CAK/Sobolev point-star line-force (O in a smooth flow. 




FIGURE 3. Left: Results of ID Smooth-Source-Function (SSF) simulation of the line-deshadowing 
instability. The line plots show the spatial variation of velocity (upper) and density (lower) at a fixed, 
arbitrary time snapshot. The corresponding grey scales show both the time (vertical axis) and height 
(horizontal axis) evolution. The dashed curve shows the corresponding smooth, steady CAK model. Right: 
For 2DHH-1DR SSF simulation, grayscale representation for the density variations rendered as a time 
sequence of 2-D wedges of the simulation model azimuthal range = 12° stacked clockwise from the 
vertical in intervals of 4000 sec from the CAK initial condition. 
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FIGURE 4. Greyscale rendition of the evolution of wind density and temperature, for time-dependent 
wind-instability models with structure formation triggered by photospheric perturbations. The boxed 
crosses identify localized region of clump-clump collision that lead to the hot, dense gas needed for a 
substantial level of soft X-rays emission. 



material. As noted for the pure- absorption models, this anti-correlation between velocity 
and density arises because the unstable linear waves that lead to the structure have an 
inward propagation relative to the mean flow. For most of the wind mass, the dominant 
overall effect of the instability is to concentrate material into dense clumps, which can 
lead to substantial (factor 3-5) overestimates in the mass loss rate from diagnostics (e.g. 
radio or Balmer emission) that scale with the square of the density. 

The presence of multiple, embedded strong shocks suggests a potential source for the 
soft X-ray emission observed from massive star winds; but the rarefied nature of the 
high-speed gas implies that this self-excited structure actually feeds very little material 



through the strong shocks needed to heat gas to X-ray emitting temperatures. To increase 
the level of X-ray emission, Feldmeier [8] introduced intrinsic perturbations at the 
wind base, assuming the underlying stellar photosphere has a turbulent spectrum of 
compressible sound waves characterized by abrupt phase shifts in velocity and density. 
These abrupt shifts seed wind variations that, when amplified by the line-deshadowing 
instabilty, now include substantial velocity variations among the dense clumps. As 
illustrated in figure HI when these dense clumps collide, they induce regions of relatively 
dense, hot gas which produce localized bursts of X-ray emission. Averaged over time, 
these localized regions can collectively yield X-ray emission with a brightness and 
spectrum that is comparable to what is typically observed from such hot stars. 

Because of the computational expense of carrying out nonlocal optical depth integra- 
tions at each time step, such SSF instability simulations have generally been limited to 
just ID. More realistically, various kinds of thin-shell instabilites [|39l] can be expected 
to break up the structure into a complex, multidimensional form. A first step to mod- 
elling both radial and lateral structure fsl] is to use a restricted "2D-H-I-1D-R" approach, 
extending the hydrodynamical model to 2D in radius and azimuth, but still keeping the 
ID-SSF radial integration for the inward/outward optical depth within each azimuthal 
zone. The right panel of figure [3] shows the resulting 2D density structure within a narrow 
(12°) wedge, with the time evolution rendered clockwise at fixed time intervals of 4000 
sec starting from the CAK initial condition at the top. The line-deshadowing instability 
is first manifest as strong radial velocity variations and associated density compressions 
that initially extend nearly coherently across the full azimuthal range of the computa- 
tional wedge. 

But as these initial "shell" structures are accelerated outward, they become progres- 
sively disrupted by Rayleigh-Taylor or thin- shell instabilities that operate in azimuth 
down to the grid scale = 0.2°. Such a 2DR-I-1DH approach may well exaggerate 
the level of variation on small lateral scales. The lack of lateral integration needed to 
compute an azimuthal component of the diffuse line-force means that the model ignores 
a potentially strong net lateral line-drag that should strongly damp azimuthal velocity 
perturbations on scales below the lateral Sobolev length Zq = rvth/vr [36]. Presuming 
that this would inhibit development of lateral instability at such scales, then any lateral 
breakup would be limited to a minimum lateral angular scale of A0m,„ ^ Iq/t = Vth/vr ^ 
0.01 rad ~ 0.5°. Further work is needed to address this issue through explicit incorpora- 
tion of the lateral line-force and the associated line-drag effect. 

Overall, both the ID and 2D SSF simulations suggest that the LDI results naturally 
in strong, small-scale clumping, implying a significant (factor several) reduction in 
mass loss rates deduced from density- squared diagnostics. On the other hand, the small 
scale makes it unlikely individual clumps can become optically thick, and so lead to a 
signifcant wind "porosity" that would weaken single-density processes like bound-free 
absorption of X-rays [|2lL ll7n. Together these suggest the observed weak asymmetry of 
X-ray emission lines from massive stars may be mainly the result of a reduced mass loss 
rate, and not wind porosity. 



Escape Integral Source Function Method and Wind Initiation 



Finally, a full linear analysis of the effects of scattering on the line-deshadowing 
instability [[26|] shows that perturbations in the line source function can lead to unstable 
wavemodes that now have outward phase propagation. Since such outward wavemodes 
have a positive correlation between density and velocity, this raises the possibility that 
the nonlinear amplification of velocity perturbations might now lead to strong forward 
shocks; this could cause a much larger fraction of wind material to be heated to X-ray 
emitting temperatures, as needed to explain the observed X-ray emission i32ll . 

These possibilities motivated development [23] of a new Escape Integral Source 
Function (EISF) that now re-uses the nonlocal escape probabilities b± defined in eqn. 
(fT9l ) to estimate the dynamical variation of the source function, 

,J;^;l^^y (20) 

Using this S{r) within a formal solution for the diffuse intensity thus provides the basis 
for a computation of the diffuse line-force that accounts for both the line-drag and phase 
reversal effects [23]. 

Application in an EISF simulation code [i24l] indeed confirms that the initial linear 
onset of the wind instability is dominated by variations with a positive correlation 
between density and velocity. However, the self-shadowing of radiative driving from 
such denser, faster structures makes their growth saturate at a low amplitude, about the 
thermal speed Vfh- Eventually, the wind structure thus again becomes dominated by the 
slower, but less- saturated growth of variations with the usual anti-correlation between 
velocity and density. Overall, the fully developed wind structure in EISF models thus 
turns out to be quite similar to the SSF case, with no apparent tendency to produce 
denser, forward shocks that might give stronger X-ray emission. 

But perhaps the biggest surprise from these EISF simulations regards the subtle nature 
of the wind initiation near the sonic point base. Because the ion thermal speed Vth is only 
moderately smaller than the sound speed a, the Sobolev length Igob = Vth/ {dv/dr) near 
the sonic point becomes comparable to the characterize density/velocity scale length 
// ~ v/ {dv/dr), thus raising questions about the suitability of using the Sobolev -based 
CAK line-force to model the wind driving in this transonic region. The left panel of 
figure[5]plots the spatial variation of various line-force components within the transonic 
wind based of a CAK initial condition, for both SSF (upper) and EISF (lower) models 
with a typical thermal to sound speed ratio Vth/a = 0.28. Note first that in both models, 
the direct force lies well below the CAK/Sobolev force throughout the subsonic region; 
as first noted by Poe et al. [jH]], this gives non-Sobolev pure-absorption models a quite 
different character than the usual CAK/Sobolev solution. 

However, in the SSF model with scattering, there now develops in this region a 
positive diffuse force, resulting from the strong steepening of the velocity gradient 
near the sonic point, which makes b > and so through eqn. (fTSl) gives a positive 
diffuse force. Quite remarkably when added to the direct term, this gives a total force 
that follows very closely the CAK/Sobolev force throughout the subsonic and transonic 
region. This explains why the (drag-effect stabilized) base outflow in SSF simulations 
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FIGURE 5. Left: Spatial variation of the velocity and various components of the line force near the 
transonic wind base for the SSF (top) and EISF (bottom) models, computed at the CAK initial condition. 
The location of the sonic point is indicated by dotted lines in the uppermost panel. Right: a. Radial 
variation near the transonic wind base of the velocity v and the line-ensemble-averaged source function 
Sa (normalized to the optically thin value). The dashed curve applies at the CAK initial condition, while 
the solid curve gives results for at time 100 ksec later in an EISF simulation, b. Corresponding variations 
of this scaled source function Sa vs. wind velocity. 



follows very closely the CAK solution, and why co-moving-frame models [30i\ of 
line-driven winds also match well the (finite-disk-corrected) CAK model. In effect, 
the addition of the SSF diffuse force corrects for the non-Sobolev error in the pure- 
absorption force, making the total force accurate to a higher order in the parameter 
hob/H. This suggests that the validity of a CAK/Sobolev approach, for which the diffuse 
force is normally assumed to vanish and thus play no role, actually depends crucially on 
the scattering nature of line-transport for proper wind initiation near the transonic base! 

But the situation is further complicated by the EISF results, which instead give a neg- 
ative diffuse force in the subsonic region. This stems from the strong "dip" in the EISF 
source function near and below the sonic point, as shown in the right panels of figure [51 
Much like the SSF diffuse force, this dip results from the velocity gradient steepening 
and associated asymmetries in the escape probabilities near the sonic point. The lower 
subsonic gradient implies a lower to shift core radiation into line resonance, while the 
large outward gradient leads to a higher b- that enhances the escape from the resonance 



zone. This lowers the mean intensity, and by eqn. (I20l) . the source function. But then the 
strong outward increase in mean intensity from this dip implies a strong inward dijfusion 
of line-flux, giving rise to the inward diffuse line-force. The net result is to weaken the 
total force even below its pure- absorption value, yielding then a reduction in the mass 
flux that can be driven through the sonic point. 

While the EISF method does not represent a fully consistent radiative transfer solu- 
tion, such source function dips have been seen in model atmospheres with full NLTE 
solutions of the line transfer Il37n. It is unclear why the effect was not seen in the co- 
moving-frame calculations of OOll . but one possibility is that this analysis was for a 
relatively dense supergiant wind, for which coupling to the continuum and thermaliza- 
tion effects might help keep the line source function elevated, and smooth, in the sonic 
region. If true, this would imply that the SSF method and results might be quite appropri- 
ate for dense winds from OB supergiants. However, for late O and cooler main sequence 
stay, the effects of a source-function dip could be significant, perhaps even playing a role 
in their "weak wind" problem ||33|]- Such effects should also be considered in the context 
of recent renewed interest in the role of atmospheric microturbulence on the wind mass 
loss rate [fl4l] . 

In conclusion, it seems that some 40 years after the initial work on line-driven stel- 
lar winds, both wind instability and wind initiation remain key issues for understanding 
their radiation hydrodynamics. To end with a personal note, I first learned about stellar 
winds as a graduate student in Dimitri Mihalas' Spring 1977 course on "Stellar Atmo- 
spheres" at the University of Colorado, and indeed even wrote my term paper on the 
subject. In grading that paper, Dimitri wrote on the cover page, "Interesting paper; I 
hope you will follow it up with some future research". In perusing this volume of papers 
from the conference in honor of his 70th birthday, perhaps he will be amused to see that 
at least one ex- student has endeavored to take his advice. 
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